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INTEGRATING MATRICES FOR ARBITRARILY 


SPACED GRID POINTS 

By 

william D. Lakin^ 

SUMMARY 

Integrating matrices form the basis of an efficient numerical 
procedure for solving differential equations associated with rotating 
beaun configurations. By expressing the equations of motion in matrix 
notation, utilizing the integrating matrix as an operator, and applying 
the boundary conditions, the spatial dependence is removed from the 
governing partial differential equations and the resulting ordinary 
differential equations can be cast into standard eigenvalue form. 

Previous derivations of integrating matrices are limited to the case of 
equally spaced grid points and approximation by interpolating polynomials. 
The restriction to equally spaced grid points may not be appropriate for 
some beam configurations. This report derives integrating matrices for 
arbitrarily spaced grid points using either interpolating or least- 
squeures fit orthogonal polynomials. Several previously unnoticrd 
features of the equally spaced grid case are also discussed. 


INTRODUCTION 


The equations of motion governing the vibrations and aeroelastic 
stability of such rotating structures as helicopter rotor blades and 
propeller blades have no closed form solution, and approximate methods 
of solution such as asymptotic techniques, Galerkin's method, or direct 
numerical integration must be employed. A numerical procedure based on 
the use of integrating matrices (refs. 1, 2) has been employed to solve 
for the vibrations and stedjility of a wide variety of rotating beam 
configurations (refs. 3, 4). The integrating matrix provides a means 
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for numerically Integrating a function that is expressed in terms of the 
values of the function at increments of the independent variable. By 
expressing the equatiors of motion in matrix notation, utilizing the inte- 
grating matrix as an operator, and applying the boundary conditions, the 
spatial dependence is removed from the governing partial differential 
equations and the resulting ordinary differential equations can be cast 
into standard eigenvalue form. Solutions can now be determined by 
standeurd methods. Historically, integrating matrices have been derived 
for equal increments of the independent variable by expressing the 
integrand as a polynomial in the form of Newton's forward-difference 
interpolation formula. However, many practical beam configurations have 
spanwise variations in the sectional properties which require the use of 
an integrating matrix which can treat unequal increments. Thus, a need 
exists for integrating matrices which can accommodate arbitrary increments 
in the indeperient variable. 

This report documents the derivation of integrating matrices which 
are valid for arbitrary increments of the independent variable. To set 
notation, the section titled "Existing Theory for Equally Spaced Grid 
Points" begins with a brief review of existing theory for equal incre- 
ments in the independent variable. The integrating matrix is derived 
here tnrough use of an interpolation polynomial. In addition to the 
basic derivation, this section discusses several previously unnoticed 
features of this approach. Integrating matrices for arbitrary increments 
of the independent variable are derived under "Theory for Grids with 
Arbitrary Spacing." Methods based on both interpolation and least- 
squares polynomial fitting are discussed and related. Finally, to test 
acctiracy, the last sectio.*^ of this report ("Verification of Accuracy 
Using Integrating Matrices based on Grids with Unequal Spacing") describes 
the application of the generalized integrating matrices to several problems 
whose solutions are known. A nimiber of sample matrices are given in 
Appendices A and B. 
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EXISTING THEORY FOR EQUALLY SPACED GRID POINTS 




The integrating matrix method is a technique for obtaining approxi- 
mations to the integral of a continuous function on a finite interval. 

To illustrate the basic ideas in the method, and to develop notation, 
this section will briefly review the existing theory for equally spaced 
grid points. 

Let f (x) be a continuous function on an interval [a, b] . Suppose 
the values of f (x) are known at the N+1 equally spaced grid points 

XQ = a , Xi = XQ + ih (i = 0, 1, . . . , N) (1) 

where h = (b - a)/N and x^ = b. For convenience, we also define 


f. = f(x.) (2) 

1 X 

Information on the function at the grid points may now be used to approxi- 
mate the function on each of the N subintervals [x^ , (j = 0, 

N-1) in the interval [a, b] . Consider, for example, the subinterval 
[XQ, xj]. Assume that f(x) can be approximated by an nth degree poly- 
nomial (n S N) on each subinterval and write f (x) - go (p) where go(p) 
is Newton's nth degree forward-difference interpolation polynomial 

9o(P) = fo + P(P - + 37 P(P - 1) (p - 2)A3fo 

+ . . . + -T P(P - • • • (P - n + DA^fo (3) 

n i 

In this formula 


P = 



(4) 


3 



while 




(E - 


D’^fo 


witl 


E^fo 


f-i, (j = 1, .. 


n) 


( 5 ) 


i.e. , E is a shift operator. 

As dx = hdp, f (x) = (p) implies that for the integral of f (x) 

over the interval (xq# xj] 



1 

f(x)dx = h 

0 



go (p) <ip 


( 6 ) 


Integration of equation ( 3 ) with respect to p and considerable 
manipulation now give the approximate result 


f(x)dx = h[aoofo + ^01^1 + •*• + aon^nl 

*0 

For example, if n = 3 , the approximation to the integral in equation 
( 7 ) is 

f(x)dx = ( 9 fo + 19 fi - 5f2 + fa) (8) 

Approximations to the integral of f(x) on each of the subintervals 

IXj, ^j+2^ (0 E 3 S N- 1 ) may be similarly obtained. Using the appropriate 

set of n+1 grid points, say x , x x , , which includes x. 

m m+1 m+n 3 

and f (5c) is approximated by an nth degree Newton's forward-differ- 

ence interpolation polynomial (p) . This polynomial is then integrated 
with respect to p from j to j+1 and the result manipulated into 

the form of a linear combination of the function values f , . . . , f ^ 

ni m+n 
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j+1 j+1 

J f(x)dx = h 

3 


J ' •• 

/ g. (p)dp = h[a. f + a. .,f 

j 3B1 m 3nri-l m+1 


X . 

3 


+ . . . + ®jjn+n^tn+n^ 


( 9 ) 


For example, if n » 3, approximations to the integrals of f(x) over the 
subintervals [xijXo], [x„ ,, x„] are 

N 


/■ 


f(x)dx + 13fl + 13f2 - f3) 


/• 


f(x)dx = + 13f3 - f 4 > 


Vi 


I 

V2 


f (x)dx 


> h 


~(-f„ , + 13f , + 13f , - f ) 

24 N-3 N-2 N-1 N 


I 

/ 


£(x)dx = 


“n-2 * ^’Vl " 


N-1 


(10) 


Having obtained approximations to the integral of f (x) over each 

b 

of the M subintervals in [a, b] , an approximation to I £(x)dx itself 

b J 


is now easily obtained. We need only note that J' f(x)dx can be written 

a 

as the "collapsing sum" 
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f (x)dx 


( 11 ) 


b 


I 


f(x)dx = 



+ 


X2 


/ 


XI 


+ 


+ 



and hence sum the approximations on the N subintervals. 


The integrating matrix representation for the approximation to the 
integral of f (x) uses the above information but puts it in a compact 
matrix form suitable for matrix manipulation Let the superscript T 
ci* a vector denote the treinsposei and define the N+1 dimensional 
column vectors {f} and {f} by 


{f} 




( 12 ) 


and 


(F) 




f (x)dx. 



(13) 


Then, relations (7) and (9) for all j = 0, ... , N-1 can be consolidated 
into a single matrix equation 


{f} = [A ]{f} (14) 

n 

where [A ] is an N+1 x N+1 matrix. As the first element of {f} is 0, the 
n 

first row of [A ] contains all zeros. Each of the remaining N rows of 
n 

[A ] contains a group of n+1 nonzero elements. For example, by equation 
n 

(7), the second row of [A ] is 

n 

aoo ^01 ^02 ••• ^On ^ •••• ^ 

Several sample matrices [A^] are shown in Appendix A. Finally, in accord 
with the collapsing sum [eq. (11)], the integrating matrix [I^] is 
obtained by left-multiplying both sides of relation (14) by the N+1 x N+1 
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suiaaing matrix 



i.e. ■ 1 if i i j but b^^ « 0 if i < j. If {f} is the N+1 
dimensional colvimn vector 


(15) 




f(x)dx, 



f(x)dxj 


(16) 


we now obtain the integrating matrix relation 


{P} - 


(17) 


with 


[I 1 » [BH?. 1 (18) 

n n 

As [B] is known a priori, the derivation of the integrating matrix 

(I ] is equivalent to the derivation of the matrix [A ] in equation (14). 
n n 

We also note that while [A ] depends on the number euid spacing of the 

n 

grid points in the interval of interest and on the degree of the approxi- 
mating polynomials, it does not depend on the values that the function 

tadces at the grid points. This separation of grid dependence and function 
dependence is one of the major strengths of the Integrating matrix. 

As described above, the theory developed in reference 2 is based on 
the use of a Newton's forward-difference interpolation polynomial to 
approximate the function on each subinterval. However, as the dimension- 
less vari<d>le p is based on the unlfona spacing h, it is impossible 
to generalize this derivation directly to grid points with unequal 
spacing . 
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Tvk> subtle aspects of the above derivation of the integrating 
matrix for equal spacing have apparently not been recognized. One 
point involves the use of Newton's forward-difference interpolation 
formula rather than a Lagrange interpolation formula. The former 
formula is equivalent to consolidating powers of x with coefficients 
that involve combinations of the functional values, i.e. 

f(x) ■ oq + oix + ... + (19) 

where uq, ... , are linear combinations of the f^. A leirge 
amount of manipulation may be required to turn equation (19) into an 
explicit linear combination of the functional values, which is the 
form that is really required to derive the integrating matrix. By 
contrast, the Lagrange interpolation formula is already in tlie proper 
form without manipulation, i.e. functional values are isolated and 
coefficients involve linear combinations of powers of x. For example, 
the Lagrange interpolation polynomial of degree n on the n+1 points 
K-), x^ , ... x^ is 

P(x) * foito(x) + fiHi(x) + ... + fnfn<x) (20) 

where 

(x-Xq) ... (X-Xj_l) (X-Xj+l) ... (X-Xjj) 

^3^*^ ~ (Xj-xo) ... (Xj-Xj. j^) (Xj-Xj^j^) . . . (xj-x„) 

It is certainly somewhat easier to integrate equation (19) than equations 
(20) and (21); however, equation (20) is clearly in the precise form re- 
quired by the integrating matrix derivation. Further, unlike Newton's 
difference interpolation formula, Lagrange interpolation is not intrinsi- 
cally limited to equally spaced grid points. 

A second important point involves the choice of a forward-difference 
formula. Rather than using forward differences, which involve starting 
at xq and moving to the right across the interval to x^, one oould 
equally well use backward differences, start at x„, and move to the 
left toward xq. The two choices are equivalent when the number of 
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interpolation points n-t-1 is even as, away from the end points, the 

interval [x., x.^.J is centered in the interpolation interval [x^, 

3 3+1 Y Y+n 

where y * 3 ■ (n-l)/2 (see fig. la). When n+1 is odd, however, two 
distinct choices axe possible (see figs, lb and Ic) . Use or forward 
differences leads to a "left biasing" with one more interpolation grid 
point to the right of left of x^ (fig. lb) . 

Similarly, backward differences lead to a "right biasing" (fig. Ic) . As 
the sample matrices in Appendix A for n « 4 and n ■ 6 indicate, the right 
and left biased matrices are related, but quite different. In particu- 
leur, a different row "marches" across the center of the matrix [A^] . 

Hence, even for equally spaced grid points, care must be taken when an 
odd number of interpolation grid points is used to approximate the 
function on each subinterval. No general rule can be formulated for 
preferring left to right biasing when n+1 is odd. The choice of 
which . ype of biasing, if any, to use must be made in the context 
of each 'rxdividual application. 

THEORY FOR GRIDS WITH ARBITRARY SPACING 

Let XQ, Xj, ... , Xjj be N+1 grid points such that 

a » x ,3 < xi < ... < x^ ■ b (22) 

Also, let f(x) again be a continuous function on the interval [a, b] 
whose values f^ « f(x^) are known at the grid points. In this section, 
we do not assume that the grid points are equally spaced in [a, b] . 

The key step in the derivation of the integrating matrix is use of 
the data at the grid points to obtain an appropriate approximation to 
the function f(x) in each subinterval [x^, (3 “ 0» ••• » N-1) . 

For arbitrary grid spacing, these approximations, and the corresponding 
integrating matrices, can be obtained in several ways. Two methods, 
interpolation and l.’iast-squares polynomial fitting, are discussed in thi's 
report. 
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(a) No biasing 



(b) Left biasing 



(c) Right biasing 


Figure 1. Examples of grid points used to approximate the 

function on the interval Ix^, Approximating 

points are marked with O. In (a), n+1 « 4, 
while in (b) and (c) n+1 * 5. 
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Interpolation 


Both the classical Lagrange interpolation formula and the more 
modem divided difference interpolation formulas are valid for arbi- 
trarily spaced grid points. Indeed, the Newton's difference formulas 
are the specialization of divided differences to uniform grid spacings. 
Unfortunately, like the Newton's difference formulas, divided differences 
in effect consolidate powers of x. The resulting coefficients again 
involve linear combinations of the data and are not well suited to use 
in deriving the integrating matrix. For thi' reason, we shall restrict 
attention here to the more appropriate Lagrange interpolation formula. 


The first step in approximating f(x) by an nth degree interpo- 
lating polynomial on the interval [x^, » j “ 0» ••• » N-l» is to 
determine the appropriate set of n+1 grid points 


G 

y 



Vi' 



(23) 


to use in the Lagrange interpolation formula. The index y in equation 
(23) will be a function of both j and n, and will fxurther depend on 
whether n+1 is even, odd with right biasing required, or odd with 
left biasing required. For grids xq/ ••• , with n-t-1 points, there 
will be N - n + 1 different sets of n-t-l consecutive grid points. 

Define Yj(j, n) , Y2(j» Y 3 (j» n) by the relations 

Yi (j. r ) ■ j - (24) 

Y2(j» n) - j - I (25) 


and 


Y3(j» n) ■ j - f + 1 


(26) 


Then, with 
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if n+1 

is even 


Y(j» n) ■ 

Y2 

if n+1 
biasing 

is odd with right 

(27) 


Y3 

if n+1 
biasing 

is odd with left 



the appropriate ;!^8finition of the interpolation set index Y(jrn) in 
equation U3) is 


0 

if 

Y 

< 

0 


Y(jr n) * { Y 

if 

0 

< 

Y < N-n 

(28) 

1 N-n 

if 

Y 


N-n 



Having properly defined the interpolation index y, and hence the 
interpolation set the nth degree Lagrange interpolation polynomial 

for f(x) on the subinterval (x^ , ' j ® 0, , N-1, is of the 

form 


f(x) = f (x) + f + ... + f ^ 

Y Y Y+1 Y**"! Y+n Y"*"^ 


(29) 


where the Lagrange coefficients £^"^Nx) are now 

m 


, , (x - X ) . . . (x - X , ) (x - X ^ , ) . . . (x - X ) 

I (Y) ^ Y mil m+1 1±£L__ 

m (x - X ) . . . (x - X , ) (x - X ^,) ... (x - X ) 


(30) 


m Y 


m m-1 m m+l 


m y^n 


If we further define 


F 


(j) . 


m 



(x)dx 

m 


(31) 


then equation (30) immediately shows that the (j+2) row of tho matrix 

[A ] for j ® 0, ... , N-1 is 
n 
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0 ... 0 ... 0 ... 0 (32) 

_ y y*l Y+n , 

Y zeros N-n-Y zeros 

As before, the first row ‘.f [A ] contains all zero elements. 

n 

The derivation of the integrating matrix via Lagrange interpolation 

is quits straightforward. However, in practice, the actual computation 

of the nN-fN quantities in the matrix [A ] for arbitrary grid 

n n 

spacings may be algebraically difficult. This is because the Lagrange 
(y) 

coefficients i (x) in equations (30) and (31) involve a product of 

IR 

n factors rather than simply powers of x. Further, evaluation of the 

(y) 

denominators in the £ (x) requires confutation of a large number of 

xn 

differences of the form k and | k - m | s n) . Fortunately, 

these implementation problems are not associated with the next method to 
be considered. Further, the following method for deriving an integrating 
matrix includes interpolation as a special case. 

Least^Squares Polynomial Fitting 

In interpolation, the interpolating polynomial agrees exactly with 
the function at the n+l interpolation grid points. However, in certain 
situations, interpolation may not make the best use of the data at the 
grid points. Suppose, for exanfle, that the degree of reliability of 
the values assigned to the function at the grid points is not well esta- 
blished. For example, this might be the case if the data on the function 
comes from experimental measiurements and small random measurement errors 
are suspected. In these circumstances, as pointed out in reference 5, 

" ... it is foolish (and, indeed, inherently dangerous) to attempt to 
determine a polyn<xnial of high degree which fits the vagaries of such 
data exactly and hence, in all probability, is represented by a ciurve 
which oscillates violently about the curve which represents the true 
function." Fitting the data with a lomr degree polynomial may be far 
more appropriate as the resulting approximation to the function will 
"smooth" out unrealistic oscillations. 
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In the present context, we will consider approximating f (x) on 
each subinterval by a least- sgueures polynomial fit based on data at 
n+1 grid points with polynomials of degree k (where k S n) . To briefly 
review least-square polynomial fits, let po (x) , Pi (x) , ... , 
a sequence of k-t-1 appropriately chosen polynomials of degrees 0, 1, 

2, ... , k. Let bhe kth degree polynomial 

k 

m=0 

Suppose we wish to approximate f (x) on the interval [ao# aj^] 
containing the n+1 grid points ao^ ai, a2# ••• # Then (x) 

is the required least-squares polynomial fit of degree k on this 
grid if the sum of the errors squared over the grid is a minimum, i.e. 

n k 2 

£=0 ^ m=0 

We note that if k = n, the interpolating polynomial on the grid aQ, aj, 
... , a^ has the minimvim value of e as it agrees with f (x) exactly 
at all grid points. Hence, interpolation can be considered a special 
case of least-square polynomial for which k = n. 

The coefficients 6 in the kth degree least-squares fit f, (x) 
m K 

are solutions of the k+1 linear ecjuations =0 (m = 0, ... , 

k) . Solution of these equations is considerably simplified if the 
polyncxnials in chosen to be orthogonal on the grid 

points ag, ... , a^. In the present context, this means 


n 

/ ^ p (a^^) P^<aj) *0 for m s 

£*0 


(35) 
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Orthogonal polynanials on a discrete grid nay be defined recursively. 
Let the orthogonal polynanials on the grid be denoted by (x) and 
let 


n r -|2 

" 2 *i[^j ^®i^J 

i=*0 

(36) 


(37) 

y... * T./S. and v. = S /S. , 

D+1 7 1 3 j 3-1 

(38) 

Po(x) = 1 

(39) 

Pl(x) = X - yi 

(40) 


and 


Pj^l(x) * (x - 

The coefficients 0 in the least-squares orthogonal polynomial fit of 
xn 

degree k s n on the n+1 grid points ao r ... , a,^ are now 



(42) 


where y^ * f (a^) . 
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Equations (33) and (42) as written involve linear combinations of 

the grid data in each coefficient B and so must be modified for the 

m 

present purposes. However, the required modifications to isolate grid 
data are quite straightforward. We obtain the basic relation 

n 

i=0 

where the coefficients q^(x) are given by 



Equations (43) and (44) may now be used to approximate the integral of 
f(x) on any interval contained in [ao, ^nl • note that the only 
integrals contained in such approximations will be integrals of the 
orthogonal polynomials (x) (j = 0, ... , k) which are quite easy to 
evaluate. Hence, the implementation difficulties associated with 
Lagrange interpolation will not be present when the least-squares fit 
is used. 

The above discussion gives the least-squares orthogonal polynomial 
fit to f (x) based on data at the grid points aof aj, ... , aj^. We 
now wish to apply these results to approximate the integral of f (x) on 
the subintervals [x^, (j = 0# ••• r N-1) . For this purpose, let 

the set be as defined in equations (23) through (28) and let 

ao — Xyr aj = Xy+lr • • • 9 ^n ~ ^Y+1 


We also denote by (x) (m = 

m 

of degree less than or equal to 


0, ... , k) the orthogonal polynomials 
k on the grid G^. As the grid data 


( 43 ) 


(44) 


( 45 ) 
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is now, by equation (45), y^ » by equations (43) 

and (44) the appropriate fit for f(x) on the interval [x., x..,] is 

3 3+1 

now 


f (X) = (X) 




i«0 


where 


( 46 ) 


^ ^*Y+1^ 

(y)(x) = > -2 I±i E 


qi'YMx) 

m*0 


£=0 ^ 


(47) 


Thus, if 




Pj+l k 

= J q {^^ x ) dx =^ 

X. hfO 


rj+1 

P^‘'^^(x..) I P^^’(x)dx 
m Y+x J m 


i : 

£*0 •* 


(48) 


equation (46) inou idiately gives for the (j+2) row of the matrix lA ] 

n 


0 

'^Y S+1 


Q. 


(j) 

Y+n 


Y zeros 


N-n-Y zeros 


(49) 


We note that, in this notation, the subscript n on the N+1 x N-tl 
matrix l A ] denotes that n-f-1 grid points are used to approximate f(x) 
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on each subinterval (x^, that the approximating ix>lynomial 

is of degree n. As noted previously, if the degree k of the 

approximating polynomial does equal n (its maximum possible value) , 

we recover the matrix [A ] obtained from interpolation. 

n 

An example of a matrix lA^J obtained via leastr-squ^res fitting is 
given in Appendix A. It is interesting to note that when k « n ^ 1, 
the middle rows of the matrices {A^] are the same as for the interpo*- 
lation case when k « n. This feature disappears when k < n ^ 1, 
Additional examples in Appendix B give matrices lA^J for various sample 
nonuniform grids. 


VERIFICATION OF ACCURACY USING INTEGRATING MATRICES 
BASED ON GRIDS WITH UNEQUAL SPACING 

As a test of accuracy, integrating matrices based on grids with 
unec[ual spacing were employed in the analysis of four types of beam 
problems previously analyzed in reference 4. These include the axial 
vibrations of a rotating beam, buckling of a rotating beam, vibration 
and stability of a rotating preconed beam, and the stability of an 
asymmetric shaft. For the present computations, the length of the 
beam was divided into 10 subintervals of unequal length (N=10) . Eight 
grid points were used to obtain approximations on each subinterval 
(n*7). Both interpolation (k=n) and least*squares polynomial fit 
(k<n) cases were considered. 

Numerical results using matrices based on seventh-degree inter- 
polation polynomials (k=n) were found to be in excellent agreement 
with results in reference 4. These earlier results were obtained 
using integrating matrices based on equally spaced grid points and 
seventh-degree polynomials. 

Most applications will involve the use of approximating polynomials 
of nMocimum degree (k«n) . However, if smoothing of experimental data is 
a factor, a least-squares polynomial fit (k<n) rather than inter-polation 
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may be «^propriate. Limited numerical evidence in the present study 
suggests that if a least-squares fit is used, the degree of the approx- 
imating polynomials should be )c>n-l. Results for k<n-l may be suspect 
if k is relatively small compared to n. 


19 


APPENDIX A 


INTEGRATING MATRICES FOR EQUALLY SPACED GRID POINTS 

This appendix contains examples of matrices [A ] for equal grid 

n 

spacings based on approximating polynomials of degree n. When n is 

odd, the resulting matrix tA^] is unique for a given grid. The matrix 

[A7] shown in table A1 is an example of this case. When n is even, 

however, two different matrices [A ] are possible for a given grid, de- 

n 

pending on whether the approximating scheme is left- or right-biased 

(see fig. 1). This fact has not been previously recognized. Tables A2 

to A5 show that the two possible matrices [A ] are related. However, 

n 

as can be seen, a different row "marches” across the center of the 
matrix in the left- and right-biased cases. It should be noted that 
the derivation in reference 2 corresponds to a left-biased approxim.-- 
tion. 

The matrix [Ay] obtained using least-squares fit polynomials of 
degree six rather than interpolating polynomials of degree seven is 
shown in table A6. 
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Table Al. The matrix {A 7 I for grids with equal spacing ~ *j “ (j=0» ••• » 10), based on seventh 

degree interpolation polynomials (k = n =« 7 ) . 
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The matrix |«„1 for ,ridx with equal upacing x^^^ - x^ = h, (J = 0 10), baaed on left 

biased approximations using fourth degree interpolating polynomials Oi - n - 4) . 
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Table A3. The matrix (A 14 ] for grids with equal spacing ••• » 10) based on right- 

biased approximations using fourth degree interpolating polynomials (k = n = 4 ). 
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Table a 5. The matrix (Agl for grids with equal spacing ~ (3=0 » ••• » 10) based on right- 

biased approximations using sixth degree interpolating polynomials (k n 6) . 
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Table A 6 . The matrix (A 7 I for grids with equal spacing = h, (j=0, ... , 10) based on least- 

squares fit polynomials of degree k - n - 1 = 6 . 
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APPENDIX B 


INTEGRATING MATRICES FOR GRIDS 
WITH ARBITRARY SPACING 

This appendix contains three examples of the matrix 
using approximation by interpolating polynomials (k^n:*?) 
unequal spacing. 


[A7] derived 
on grids with 
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Table Bl. The matrix [Ay] based on seventh degree interpolation polynomials (k = n = 7) for the grid 
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Table B3. The matrix [Ay] based on seventh degree interpolation polynomials (k = n = 7) for the grid 
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